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With sufficient data, Large Hadron Collider (LHC) experiments can constrain the selectron-smuon 
mass splitting through differences in the di-electron and di-muon edges from supersymmetry (SUSY) 
cascade decays. We study the sensitivity of the LHC to this mass splitting, which within mSUGRA 
may be constrained down to C(10~'') for 30 fb~^ of integrated luminosity. Over substantial regions 
of SUSY breaking parameter space the fractional edge splitting can be significantly enhanced over 
the fractional mass splitting. Within models where the selectron and smuon are constrained to be 
universal at a high scale, edge splittings up to a few percent may be induced by renormalisation group 
effects and may be significantly discriminated from zero. The edge splitting provides important 
information about high-scale SUSY breaking terms and should be included in any fit of LHC data 
to high-scale models. 
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TeV-scale supersymmetry is one of the most promis- 
ing solutions to the weak hierarchy problem of the Stan- 
dard Model. If TeV-scale supersymmetry is realized in 
nature, the LHC will be able to discover it by the pro- 
duction of sparticles and measure some of the param- 
eters of the minimal supersymmetric standard model 
(MSSM) . Through the use of renormalisation group equa- 
tions (RGEs) to connect the low and high scale La- 
grangians, the measured MSSM parameters can be used 
to test hypotheses, such as flavor universality (or non- 
universality) of the high-scale soft SUSY breaking terms. 
Such a link would be a welcome additional empirical win- 
dow into the unsolved problem of flavor and how it relates 
to SUSY breaking. 

Potentially one of the most accurate LHC measure- 
ments of SUSY is that of the di-lepton edge arising from 
neutralino-slepton-ncutralino SUSY cascade decays. 



If the mass ordering m^o > mj > m^o is present, the di- 
lepton mass spectrum has a prominent kinematic edge 
at 
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rrinl S> rUn, the finiteness of 



the muon mass induces fractional corrections that are 
^ {^Ibl^lusy) (where rugusy denotes the order of mag- 
nitude of sparticle masses, > 100 GeV) and may be ne- 
glected, as they are throughout this paper. Such an edge 
may be measured at LHC experiments with per-mille pre- 
cision [1]. Although the original authors commented that 
one could measure the di-lepton and di-muon endpoints 
separately [2|, it is generally assumed, motivated by the 
MSSM flavor problem, that mg 
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The flavor problem is one of the most pressing 
questions of supersymmetric phenomenology. The La- 
grangian of a TeV-scale MSSM is highly constrained by 
the absence of new contributions to low energy flavor 
changing neutral currents (FCNCs). Classic examples 
are K^K^ mixing or in the leptonic sector the fi —> 
branching ratio. These strongly suggest that new physics 
present at the TeV scale should obey to high accuracy the 
principle of MFV (Minimal Flavor Violation), at least in 
the first two generations. The /i — )■ 67 branching ratio, 
which is currently bounded at BR{ijl ej) < 1.2 x 10"^""^ 
H and in the future will be measured to at least 10"^"^ 
[4|, constrains off-diagonal propagator mixing between 
the smuon and selectron LL flavor eigenstates to 0], 

Si2,LL < 6 X 10-^ 

The RR constraints are similar over most of parameter 
space, but potential cancellations exist and so 5i2,_Ri? has 
a weaker bound than Si2,ll- 

Nonetheless, exact low-scale flavor universality is not 
expected; some non-universality will be automatically in- 
duced by RGEs and in the context of string theory mod- 
els of supersymmetry breaking loop effects are expected 
to violate flavor universality even if leading order physics 
preserves it 6]. It is therefore important to investigate 
the sensitivity of the LHC to the mass splitting between 
the right-handed selectron and smuon: 
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edges occurring at identical values. 



Am? is not directly constrained by the BR{ii 67) 
measurement as it does not give lepton flavor violation 
(LFV). LFV effects at the LHC have been studied in 
various papers 0, H, H, [lOj [HI El) [ill but are not our 
focus here. 

We may relate Am^ to the SUSY breaking terms de- 
fined at a high scale Mx through the MSSM RGEs. At 
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one-loop order [if 
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where fi is the ZJi? renormalisation scale, and m^^, 
TTip^j^), rriej^, ^-hcI rn_f/i are the soft SUSY break- 

ing masses for right- and left-handed smuons and se- 
lectrons and the Higgs field, evaluated at a renormal- 
isation scale /X. /le, are the trilinear selectron and 
smuon soft SUSY breaking couplings evaluated at a 
renormalisation scale /i. At M^, we have the bound- 
ary condition Y^{Mz) = ^/2m ^{M z) / {v cos fi) , where 
V = y/vliMz) +v^{Mz), wi,2 being the MSSM Higgs 
vacuum expectation values. We may solve the RGE for 
Y^, Y^iMx) = Y^{Mz)+0{ln{Mx/Mz)/167r^). We ne- 
gleet the ratio of the electron Yukawa coupling to the 
muon Yukawa coupling (Ye/y^) and solve eq. ([3]) to first 
order in ln(Afx/Afz)/(47r)^ to obtain 



Am^ (Mz) = Am^{Mx) 
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+K,{Mx)+mij^iMx) + 



AliMx) 



tan^/51n( ^ 
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where tan (3 — V2 jvi . We have allowed for the presence of 
a primordial mass splitting Aw? (Mx ) and used the large 
tan/3 limit in eq. (|4]) so that cos/3 « 1/ tan/9. We have 
also neglected the QED running of the muon mass below 
the electroweak scale. We see from eq. ^ that the mag- 
nitude of An? is always enhanced as tan^ /31n(Mx/Mz). 
The evaluation of Arr?(mz) is complicated by the un- 
known boundary conditions on to , w , and 
at Mx^ but within constrained models other edge mea- 
surements are expected to help pin them down. For ex- 
ample, in the mSUGRA model, the SUSY breaking scalar 
masses have a universal value of toq at Mx and the tri- 
linear SUSY breaking scalar couplings have a universal 
value of Aq. In this case, Ato^(Mx) = and the quan- 
tity in square brackets is simply 3mQ 4- A§ , which may be 
bounded from other measurements. 

From eq. [T] the variation of the edge position with the 
slepton mass is given by 
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with a fractional shift in the invariant mass edge of 
Amu Ami ( 



m^om^a 
xi xi 
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FIG. 1: The enhancement factor {Amn/mu) / {Ami/mi) as a 



function of Ami/Am2 
different values of n 



(mi — m^o)/(m^o — "i^o) for three 

X2 ' Xi 



up to terms of O ^^^^ 



We define the enhancement factor by ( " )/(^^). 

^ ^ mil ^ ' 

Fig. [T] shows the enhancement factor as a function of 
the slepton mass. As well as the mSUGRA bino- 
dominated lightest supersymmetric particle (LSP) case 
n = m^o^jmyO = 2 we also allow for a more compressed 
spectrum with different values of n. There are several no- 
table points. When mi = ^m^om-^o the shift in the edge 

vanishes to leading order in . Conversely, large split- 
tings of the di-lepton and di-muon edges can be achieved 
for relatively small splittings of the slepton and smuon 
masses. The enhancement factor is larger for more degen- 
erate masses between sparticles in the chain and can eas- 
ily be 0(10) depending upon the value of m^o /m^o . The 
enhancement diverges as the slepton mass approacWs ei- 
ther neutralino mass. The benefits of the enhancement 
may be diluted for highly degenerate spectra by the fact 
that leptons coming from such chains will tend to be 
softer and thus harder to identify and measure experi- 
mentally, and also by any change in precision due to rate 
changes from differing branching ratios. 

The muon energy calibration will depend heavily on 
their energies, but we may expect particularly accurate 
calibration when the mu edge is close to, but not on, 
the Z° pole. Electrons and muons from the Z'^ decays 
may then be used to calibrate energies and efficiencies, 
performing a small extrapolation to the relevant energy. 
For instance, the scale of the mu endpoint is 80, 118, 
99, 122, 343 GeV for the Snowmass benchmark points 
SPSla, 3, 5, 8 and 9 respectively [3. Thus 4 of the 5 
relevant endpoints are rather close in invariant mass to 
the ZO. 

To estimate the experimental sensitivity to differences 
in the positions of the di-muon and di-electron edge, 
we generate 400,000 i?-parity conserving supersymmet- 
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FIG. 2: Log likelihood fit to di-electron edge in AcerDET as 
described in the text, rriee has been placed in 0.25 GeV bins 
on abscissa. The ordinate shows events per bin per 16 fb~^. 



ric events at SUGRA Point 5 from proton-proton colli- 
sions at a 14 TeV centre of mass energy. We use HERWIG 
6.510 which reports the two-to-two (Standard Model to 
Supersymmetric Particle) cross section to be ~ 24 pb. 
Our sample of 400,000 events therefore corresponds to 
an integrated luminosity of ~ 16 ,fb^^. We pass the gen- 
erated events through AcerDET 14] which represents a 
fast simulation of a generic general-purpose LHC-type 
experiment. The main feature that AcerDET provides 
for our purposes is a reasonable model of the electron 
and muon momentum resolution for a typical LHC-type 
detector. AcerDET accomplishes this by 'smearing' the 
generated momenta by appropriate amounts dependent 
on pt, Tj and </>. AcerDET provides only minimal simula- 
tion of reconstruction efficiencies. For example, electrons 
will alway be reconstructed unless they are out of accep- 
tance or are too close to other particles. Neither does 
AcerDET model any uncertainty in the absolute energy 
scale for either electrons or muons. For this we use an 
estimate of 0.1% [ll, [ll|. Absolute energy scale calibra- 
tion at this level might appear to be a tall order - however 
the standard candles used in this calibration (Z-bosons) 
will be produced in such large numbers at the LHC that 
the estimate may even turn out to be conservative. To 
evaluate the sensitivity with which the endpoints of the 
di-electron and di-muon invariant mass spectra can be 
measured, we generate their distributions, fit them in 
the vicinity of the endpoint, and report the endpoint fit 
uncertainties. To select our event we require two oppo- 
site sign same family (OSSF) isolated leptons of greater 
than 10 GeV and missing energy greater than 100 GeV. 
At this point, the di-lepton invariant mass distribution 
should in principle end at 97.48 GeV. We fit over the 
range 60 GeV to 140 GeV. For our (log likelihood) fit we 
use a function which takes the form of a triangle distri- 



FIG. 3: Log likelihood fit to di-muon edge in AcerDET as 
described in the text, m^^ has been placed in 0.25 GeV bins 
on abscissa. The ordinate shows events per bin per 16 fb~^. 
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(GeV) 


16.0 


22145 


97.47 ±0.09 


97.56 ±0.18 


8.0 


11131 


97.41 ±0.13 


97.83 ± 0.23 


4.0 


5520 


97.54 ±0.19 


97.63 ± 0.35 


2.0 


2707 


97.52 ±0.28 


97.56 ±0.50 



TABLE I: Results of the endpoint fits for various integrated 
luminosities. The number of OSSF di-lepton events passing 
cuts and having di-lepton invariant mass below 100 GeV is 
also recorded. 

bution (i.e. f{x]e) oc xQ{e — x)Q{x) where Q{x) is the 
Heaviside step-function and e is the notional endpoint) 
sitting on top of a constant background, and with the 
entire distribution convolved with a Gaussian resolution 
of width a. In the fit, the three free parameters were 
(1) the endpoint position e, (2) its resolution a, and (3) 
the ratio of the number of signal and SUSY background 
events. The overall normalisation was fixed analytically 
to allow a log-likelihood fit to be performed with nar- 
row bins (0.25 GeV). Example fits are shown in figures [2] 
and [31 Note that we neither generate nor fit Standard 
Model (SM) events. In reality, irreducible sources of SM 
di-leptons from Z-bosons would constitute an important 
background if the di-lepton endpoint were to occur on the 
Z resonance. We ignore this special case for the moment 
because it is our primary aim to make statements about 
the "generic" di-lepton edge sensitivity. Away from the 
Z-peak, the missing energy cut will reduce the SM back- 
grounds. The results of the endpoint fits are summarised 
in Table HI 

Accordingly we parametrise the LHC di-lepton edge 
sensitivity E (which we define to be the expected frac- 



tional uncertainty in the di-muon di-electron edge differ- 
ence as estimated by the fit error on the less well deter- 
mined endpoint added in quadrature to the energy scale 
uncertainty) by 



E = \J (0.002 v/22100/7V)2 



0.0012 



(6) 



for N expected signal events in the di-lepton channel. 
In this expression, the "0.001 term" represents the 0.1% 
absolute energy scale error described earlier. Note that 
the di-lepton edge sensitivity E is not to be confused with 
the slepton mass sensitivity E defined later. 

By the above definition, the edge sensitivity S is a 
measure of the scale down to which fractional differences 
in endpoint positions Amu/mu can be measured. More 
precisely, assuming that the endpoint fit error is approx- 
imately Gaussian distributed, it should be possible to 
make an "S'l-sigma" discovery of selectron-smuon pole 
mass non-universality (i.e. rule out a null hypothesis of 
"no splitting in the di-electron and di-muon endpoints" 
at the S'l-sigma level) for a real endpoint splitting of size 
Amu according to 



Si 



Amu 



(7) 



In what follows we will therefore refer to Si as the 
"discovery significance" for selectron-smuon mass non- 
universality. We note that, when calculating Si, the 
numbers of events N^e contributing to the e+e~ di-lepton 
signal may differ from the number of /i^/i^ pairs iV^^ 
due to phase-space differences induced by the mass dif- 
ferences. If systematic uncertainties on trigger and recon- 
struction efficiencies can be controlled, this could provide 
an additional means of testing selectron-smuon mass uni- 
versality by looking at significant differences from zero in 
the statistic 



So = 



N,, - N,, 



N 



(8) 



which, like Si, will be approximately normally dis- 
tributed. We do not use 5*2 ourselves. 

In the examples which follow we calculate the TV in 
Si for 30 fb^^ of integrated luminosity at the LHC us- 
ing WIGISASUGRA1.200 and HERWIG6.5 [20]. Perform- 
ing a Markov Chain Monte Carlo maximisation on Si 
in mSUGRA, we find a maximum value 5*1 = 0.52 af- 
ter direct search constraints have been applied. Thus 
we find that the smuon-selectron splitting cannot be dis- 
criminated from zero in niSUGRA at the LHC. On the 
other hand, any significant measured difference in the end 
points at the LHC will discriminate against mSUGRA. 
A future international linear collider would achieve much 
improved accuracy [l^ upon the mass splitting and could 
be combined with other constraints to help bound tan /3 
assuming mSUGRA. 

We wish to emphasise that the reason mSUGRA fails 
to generate an observable edge splitting is not the one 
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FIG. 4: Expected 30fb~^ 1-sigma sensitivity, E, to selectron- 
smuon mass splitting in perturbed mSUGRA around SPSla. 
The region to the right hand side of the almost-vertical line 
has mf 1 < m^o . The region underneath the mostly horizontal 
line has m^o — mf < 10 GeV. The lighter lines show contours 
of logj^Q E = —2, —2.5, —3 (top to bottom) 



often suggested. It is not true that the muon and elec- 
tron Yukawa couplings are too small to play any role in 
the RGEs. Indeed, at large tan/3 the RGEs, combined 
with the enhancement factor, can generate slepton spec- 
tra giving edge splittings at the per cent level. The real 
reason is that in this case ffl is driven light and it dom- 
inates the X2 decay modes, with BR{x2 — > IrI) ^ 1 
and BR{x2 ~^ t^it) ~ 1. In models where there is extra 
third family physics that would lead to the weak scale 
mass ordering mf„ > m^o > mr > m^o, the selectron 
and smuon RGEs can be sufficient to generate electron- 
muon edge splittings that can be significantly discrimi- 
nated from zero. As an example, we consider the point 
Too = 148 GeV, M1/2 = 250 GeV, Aq = -600 GeV, 
tan/3 = 40 with mf^^{Mx) = 950 GeV. At this point, 
Amj/mj = 2.3 x 10"'^ and Amu/mu = 1.5% whereas 
S = 0.27%, allowing an (5*1 > 5)-sigma discovery signifi- 
cance for smuon-selectron pole mass non-universality. 

We now focus directly on the sensitivity to slepton 
mass splittings and analyse how degenerate the selectron 
and smuon masses can be while still allowing a 1-sigma 
sensitivity to a non-zero mass difference. We do this as- 
suming 30fb~^ in perturbed mSUGRA around SPSla. In 
perturbed mSUGRA, we take mSUGRA boundary condi- 
tions but we allow to^^ to float away from the mSUGRA 
prediction, as could be derived from Ato?(Mx) 7^ in 
Eq. m By using Eqs. [H and [71 setting Si to 1 we ob- 
tain the fractional slepton mass splitting which might be 
discriminated from zero at the 1-sigma level: 



E : 



Si=l 



{m\ - m?){mj - mK 

2 2 4 

m „m „ — mr 

Xi X2 I 



valid in the limit Amj/mj <^ 1. Fig. |3] displays E for 
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FIG. 5: Probability distribution difference in GeV"^ of SUSY 
cascade chain electrons and muons as a function of their in- 
variant mass mil (in GeV) for an endpoint of 70 GeV and a 
relative splitting of 2%. The solid line shows the distribution 
without any energy resolution taken into account, whereas 
the dashed line displays the effect of smearing due to energy 
resolution. The dotted line shows a smeared distribution with 
no mass splitting. 



a scan of perturbed mSUGRA around SPSla. The ex- 
pected sensitivity at SPSla itself is = 2.8 x 10"^. The 
strict mSUGRA prediction for the smuon-selectron mass 
splitting is lS.mi/mj — 5.9 x 10~^. We see that sensitivi- 
ties down to 0(10"'') are possible while restricting to the 
region where m^o — mj > 10 GeV to ensure sufficiently 
hard leptons. 

It is tempting to display the difference in ijl^ jJT and 
e+e" probabiUty distributions (AP) in order to look for 
a spike, see the solid curve in Fig. O In practice, en- 
ergy resolution effects would smear out this curve. In the 



dashed curve, 



and ' 



have been Gaussian smeared 



by assumed fractional resolutions of 1% and 3% respec- 
tively. We have assumed Nf.^ — N^^ and simulated no 
backgrounds. We have also plotted AP for equal end- 
points while taking the energy resolution into accovmt in 



the dotted line. Simply having a worse energy resolution 
for muons still leads to a similar feature in the difference 
near the end-point and could mislead us into thinking 
there is a splitting when in fact there is none. We con- 
clude that, in practice, the best approach will be to fit the 
di-electron and di-muon endpoints separately, and then 
examine the difference, as was done in Figs. [5] and [31 We 
note here that one could still employ the technique of 
subtracting opposite sign different flavor di-leptons from 
the di-muon or di-electron decay chain samples in order 
to subtract backgrounds (e.g. from top pairs or W pairs), 
as was recommended in the case [l| where the two sam- 
ples are summed. 

In this paper we have studied the sensitivity of the 
LHC to the jl-e mass splitting through endpoint differ- 
ences in SUSY cascade decays. Enhancement factors 
mean that measured differences in the endpoints can be 
a factor of ten more sensitive to tje mass splitting. In 
the large tan/3 limit RGE can induce edge splittings 
up to the per cent level. However, in mSUGRA the 
X2 ^ TiT branching ratio increases, significantly reduc- 
ing the X2 ^ ^R^ branching ratio. A significant rejec- 
tion of the universal smuon-selectron pole mass hypoth- 
esis would discriminate against mSUGRA. If additional 

then the RGE- 



mass terms were to cause mf 



induced edge splitting may be significantly discriminated 
from zero. The di-lepton edge splitting may be a power- 
ful measurement and discrimination tool in SUSY data 
analysis, and should not be forgotten in global fits to 
data. 
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